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Abstract. By cross-correlating the 1.4 GHz FIRST catalog of radio sources with the Abell cluster catalog we 
have found an excess surface density by a factor ~ 5 of radio sources within a projected distance, r p , of 0.1 Mpc 
(for h = 0.65) from the cluster center. The profile of the excess density can be described, for r p 3; 0.1 Mpc, by 
a /3-model with a core radius of ~ 0.70 Mpc and = 1.65. The luminosity function of cluster sources does not 
show hints of cosmological evolution over the redshift range (z S; 0.4) covered by our cluster sample. The mean 
luminosity function is in excellent agreement with the recent determination by Reddy & Yun (2004) for 7 nearby 
clusters and extends it by two orders of magnitude to higher luminosities. Its shape is very similar to that of the 
local luminosity function of field galaxies, but the space density is about 3000 times higher. When extrapolated to 
30 GHz, our luminosity function compares very favourably with an estimate obtained directly from the 30 GHz 
observations by Cooray et al. (1998). The antenna temperature contributed by radio sources within the nominal 
cluster radius of 1.7 Mpc is estimated to be ~ 13.5 fiK at 30 GHz for clusters at z ~ and decreases to ~ 3.4 fjK 
at z ~ 0.5, in the absence of cosmological evolution; it increases by a factor of ~ 1.5 within 0.25 Mpc from the 
cluster center. If the pure luminosity evolution models by Dunlop & Peacock (1990) are adopted, the radio source 
|' antenna temperature turns out to be essentially independent of redshift. 
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c/j , Key words. Cosmic microwave background — Galaxies: In this paper, we exploit the FIRST (Becker et al. 1995; 

■ clusters: general — radio continuum: galaxies White et al. 1997) survey at 1.4 GHz to investigate the 
, correlation function of radio sources with Abell clusters 

■ (Abell 1958) and to derive an estimated of the luminosity 
lj_J ! 1. Introduction function of radio sources in clusters. Coupling the FIRST 

with the GB6 (Gregory et al. 1996) survey, we obtain an 

Over the last decade the quality of observations of the estimate of the spectral index distribution that we use to 

Sunyaev-Zeldovich (SZ) effect (Sunyaev & Zeldovich 1972, extrapolate the luminosity function to centimeter wave- 

1980, 1981) towards galaxy clusters has dramatically im- lengths and to estimate the mean contamination of the 

proved, so that it is now routinely detected with high S/N SZ signal as a function of the cluster redshift. 
and even imaged (Carlstrom et al. 2002). To fully exploit 
the extraordinary information content of the effect, how- 
ever, one needs to take into account the sources of possible 
contamination. 

Since powerful radio sources are normally associated to 2.1. Source selection 
early-type galaxies and these preferentially reside in clus- 
ters, a strong over-density of radio sources, which could We have used the April 2003 update of the FIRST survey 
partly fill the SZ dip in the Rayleigh- Jeans region, is natu- catalog (sundog.stsci.edu/first/catalogs.html), containing 
rally expected in clusters. And, indeed, radio sources were 811,117 radio sources over an area of 8,422 square degrees 
found to be the major contaminant of experiments us- in the North Galactic cap and 611 square degrees in the 
ing centimeter-wavelength receivers (Cooray et al. 1998; South Galactic cap, with a detection threshold of about 1 
LaRoque et al. 2003). mJy and an angular resolution of 5". The integrated flux 
was adopted for sources labelled as extended and the peak 



2. FIRST sources within cluster solid angles 



Send offprint requests to: G. De Zotti: dezotti@pd.astro.it flux for the others. 
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We have looked for FIRST sources within a projected 
distance of 3Mpc (typical cluster virial radius), IMpc, 
and 0.25 Mpc from the Abell cluster centers, for a flat 
ACDM cosmology (fi TO = 0.3, Q A = 0.7) and H = 
6511ms- 1 Mpc" 1 (or h = H /100 = 0.65). We have kept 
only those clusters fully within the FIRST area (taking 
into account the "holes" in such area). We have also ex- 
cluded clusters partially overlapping with each other. 

Measured redshifts for 368 Abell clusters in the FIRST 
area are listed in Struble & Rood (1999); 1 additional red- 
shift measurement was found in Pimbblet et al. (2001). 
Redshift estimates for the remaining 582 clusters were ob- 
tained using the relationship between redshift and mag- 
nitude of the tenth brightest cluster galaxy derived by 
Postman et al. (1985), which includes a correction for the 
Scott effect. A comparison of estimated with measured 
redshift showed this relationship to be superior to those 
proposed by other Authors (Fullerton et al. 1972; Leir et 
al. 1974; Kalinkov et al. 1994). 

Our final sample comprises 951 clusters with a total 
of 617, 5,574, and 35,526 radio sources within 0.25 Mpc, 
1 Mpc, and 3 Mpc from the cluster centers, respectively. 
We have also selected 738 control fields, each having an 
angular radius of 0°.5 and whose centers are more than 1° 
away from the center of any of the clusters. The control 
fields contain a total of 51,770 radio sources. 




Log(r p ) 

Fig. 2. Projected two-point correlation function of radio 
sources with cluster centers, as a function of project radius 
(Mpc, for h = 0.65). The dotted line shows a /3-model fit 
with r = 0.70 Mpc and /3 = 1.65 [see Eq. ©]. 
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Fig. 1. Surface density (sr" 1 ) distribution of radio sources 
in the control fields (shaded area) and for clusters in the 
lowest redshift interval, within 0.25 Mpc from their cen- 
ters. 



2.2. Analysis of the sample 

The cluster sample has been divided into 8 redshift in- 
tervals. The mean surface densities within 0.25, 1, and 3 



Mpc from the cluster centers, and their errors are listed, 
for each interval, in Table ^ Figure ^ shows an example 
of the distribution of surface densities within 0.25 Mpc, 
no.25, compared with the distribution for control fields. 
The mean value of no. 25 (no.25) is typically almost a factor 
of 3 higher than in control fields, without any significant 
trend with redshift. Cooray et al. (1998), who surveyed 
56 clusters, within 250" from the pointing centers, at 28.5 
GHz, to a 4er detection limit ~ 0.4 mJy, found 4 to 7 times 
more sources than expected from a simple extrapolation 
of low- frequency counts outside clusters. 

The over-density towards clusters is still quite substan- 
tial within 1 Mpc from the cluster centers and is weak, but 
still statistically significant, within 3 Mpc. A Kolmogorov- 
Smirnov test gives, for the first redshift bin in Table ^ 
probabilities of ~ 2.8 x 10~ 12 and of ~ 6.5 x 10~ 8 that the 
distributions of no.25 and of n\, respectively, are drawn 
from the same parent distribution as for control fields. 
Such probability drops by orders of magnitude for the 
higher redshift bins. The difference with control fields of 
the distribution of is not significant for the lowest red- 
shift bin, is ~ 1.6 x 10~ 2 , for the second lowest redshift, is 
~ 2 x 10 -3 for the sixth bin, and is orders of magnitude 
lower for all the other bins. The densities ri\ and cor- 
relate with redshift, although the redshift dependence is 
rather weak. 

To characterize more quantitatively the over-density, 
we have computed the projected radio source-cluster two- 
point correlation function, £ p , as a function of the pro- 
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Table 1. Mean surface densities (sr _1 ) of FIRST radio sources within projected distances of 0.25, 1 and 3 Mpc 
(^0.25, n\, and 77,3, respectively) from cluster centers, for different redshift intervals (the notation xxE+y stands 
for xxx 1Q V ). N c \ is the number of clusters in a given redshift bin. The average surface density of control fields is 
(2.93 ±0.545) x 10 5 sr- x . 



log z min 


log z max 


N d 


"-0.25 


ni 


n 3 


-1.80 


-1.15 


41 


(7.7 ± 1.5)E+5 


(3.93±0.27)E+5 


(2.96±0.06)E+5 


-1.15 


-1.00 


81 


(7.8 ± l.l)E+5 


(4.15±0.24)E+5 


(3.03±0.07)E+5 


-1.00 


-0.90 


130 


(8.1±l.l)E+5 


(4.13±0.23)E+5 


(3.13±0.07)E+5 


-0.90 


-0.85 


201 


(7.7 ± 1.0)E+5 


(4.64±0.22)E+5 


(3.16±0.07)E+5 


-0.85 


-0.80 


149 


(8.1 ± 1.4)E+5 


(4.62±0.25)E+4 


(3.19±0.07)E+5 


-0.80 


-0.75 


168 


(7.8 ± 1.3)E+5 


(5.12±0.26)E+5 


(3.24±0.08)E+5 


-0.75 


-0.70 


105 


(7.9 ± 1.5)E+5 


(5.02±0.35)E+5 


(3.29 ± 0.11)E+5 


-0.70 


-0.35 


76 


(1.4 ± 0.4)E+6 


(7.10±0.61)E+5 


(3.76±0.15)E+5 



jected linear distance, r p , from the cluster center, using 
the formula: 



N c (r p ) 
Nf(r P ) ' 



(1) 



where N c is the mean number of sources within annuli 
of width Ar = 0.1 Mpc, centered in the cluster centers, 
Nf is the mean number of control field sources in areas 
equal to those of the annuli, and r p is equal to the inner 
radius of annuli plus Ar/2. The errors on £, p (r p ) where 
estimated as S£ p (r p ) ~ (1 +£ P (?>))/ (v'^Vr) [Peacock 1999, 
Eq. (16.111)]. 

As shown by Fig. [3 the surface density of sources close 
to the cluster center is about a factor of 5 higher than in 
the control fields, and decreases by more than a factor of 
10 at a radius of ~ 1 Mpc. Except for the innermost point, 
the projected correlation function, £, p (r p ) can be described 
by a /3-model (Cavaliere & Fusco-Femiano 1976) surface- 
density profile 



-(3/2)/3+l/2 



<r(r p ) 



, 1+r o / 



(2) 



with ro = 0.70 Mpc and (3 = 1.65. For comparison, Reddy 
& Yun (2004) found, for a sample of 182 radio galaxies 
in 7 nearby, rich, X-ray luminous clusters, ro = 0.74 Mpc 
(scaled to our choice for the Hubble constant) and (3 = 
1.05. 

3. Luminosity function of radio galaxies in clusters 

Reddy & Yun (2004) have estimated the mean 1.4 GHz 
luminosity function of galaxies in 7 nearby clusters con- 
taining a total of 182 NVSS (Condon et al. 1998) sources, 
assuming a uniform distribution within a cluster radius of 
1.5 Mpc for h = 0.75 (or 1.73 Mpc for the value of the 
Hubble constant, h — 0.65, adopted here). 

Since the median density of radio galaxies within 0.25 
Mpc from the cluster center is ~ 8 x 10 5 sources/sr 
(Table 0, i.e. a factor ~ 2.7 larger than the field den- 
sity, an average of 63% of the 617 galaxies within such 
radius are cluster members and we may attempt to esti- 
mate their luminosity function. Following Reddy & Yun 
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Fig. 3. Estimated luminosity functions in three represen- 
tative redshift bins: the first (stars), the third (filled tri- 
angles), and the seventh (open squares) in Tabled 



(2004), we have assumed a uniform distribution within a 
cluster radius of 1.73 Mpc (for our choice of the Hubble 
constant); the sampled volume is then the intersection of 
such sphere with a cylinder with radius of 0.25 Mpc. The 
derived densities have been scaled down by a factor of 0.63 
to take into account the foreground/background contami- 
nation and by a further factor of 4.48 to convert the mean 
density within the sampled volume into the mean density 
within the adopted cluster radius. 

To check for possible evolutionary effects we have car- 
ried out separate estimates for each redshift bin in Tabled 
Results for three representative redshift bins arc shown in 
Fig- El No indication of evolution is apparent. The low lu- 
minosity excess found for the lowest luminosity bin is likely 
due to contamination. Since we consider sources within a 
fixed physical radius the number (but not the fraction) of 
outliers (proportional to the corresponding solid angle) in- 
creases with decreasing redshift; the flux densities of field 
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Fig. 4. Present estimate of the mean 1.4 GHz luminosity 
function in clusters (open circles) compared with that by 
Reddy & Yun (2004; filled triangles) and with the local 
luminosity function in the field (Magliocchetti et al. 2002; 
crosses). 



sources are preferentially close to the survey limits and 
therefore affect more the lowest luminosity bins. 
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Fig. 5. Estimates of the luminosity function of cluster ra- 
dio galaxies at 30 GHz. The open circles show the extrap- 
olation to 30 GHz of the luminosity function in Fig. 0] 
The filled triangles represent our estimate based on the 
30 GHz data by Cooray et al. (1998). 



Table 2. Estimate of the 1.4 GHz mean luminosity func- 
tion of cluster radio sources. 





logfltWHz-'sr- 1 )] 


logKMpc^dlogL- 1 )] 




20.3 


0.12 ±0.58 




20.7 


0.43 ± 0.38 




21.1 


0.44 ±0.16 




21.5 


-0.83 ±0.12 




21.9 


-1.08 ±0.09 




22.3 


-1.14 ±0.09 




22.7 


-1.19 ±0.09 




23.1 


-1.35 ± 0.11 




23.5 


-1.60 ±0.15 




23.9 


-2.48 ±0.41 




24.3 


-2.66 ±0.50 



In the absence of evidences of cosmological evolution 
of the luminosity function, we have lumped all galaxies to- 
gether to derive the mean luminosity function presented 
in Table |3 Figure 01 shows that our estimate nicely agrees 
with that by Reddy & Yun (2004), based on totally differ- 
ent approach, and extends it by two orders of magnitude 
in luminosity. The figure also compares the cluster with 
the field local luminosity function recently determined by 
Magliocchetti et al. (2002). The shape of the luminosity 
functions is similar but density of radio sources in clusters 
is about 3000 times higher than in the field. For compar- 
ison, the matter overdensity in a cluster formed at a red- 
shift z f is ~ 200(1 + z f f ~ 3000 for a formation redshift 
z f ~ 1.5. 

The integral of the 1.4 GHz cluster luminosity func- 
tion, made adopting the estimate by Reddy and Yun 
(2004) for logL(WHz _1 sr _1 ) < 22 and our results for 
higher luminosities, gives a mean cluster radio luminosity 

log I<clustci\ 1.4GHz — 23.35. 

4. Radio source contamination of 
Sunyaev-Zeldovich signals 

To extrapolate the 1.4 GHz cluster luminosity to the 
higher frequencies used for measurements of the Sunyaev- 
Zeldovich effect, we need the spectral index distribution of 
sources. To get some information on that, we have cross- 
correlated the GB6 4.85 GHz catalog (Gregory et al. 1996) 
with our cluster sample, looking for sources within 0.25 
Mpc from the cluster center. The number of matches is rel- 
atively small (46 sources) because of the relatively bright 
flux density limit (18 mJy, to be compared with the 1 
mJy FIRST limit). The identification of GB6 with FIRST 
sources is complicated by the relatively poor GB6 angu- 
lar resolution. In fact, in several cases we have multiple 
FIRST sources within the GB6 beam (FWHM ~ 3'. 5). In 
these cases, we have taken as the 1.4 GHz flux density of 
the GB6 sources the sum of FIRST flux densities weighted 
with a Gaussian response function with FWHM ~ 3'. 5 
and axis towards the nominal GB6 position. Assuming 
power-law spectra (S cx v~ a ) we found that, of the 46 
GB6 sources, 12 are flat-spectrum (a < 0.5) and 34 are 
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Fig. 6. Mean contamination of the SZ signal (in antenna 
temperature) by radio sources as a function of cluster red- 
shift for 4 frequencies, 30, 44, 70, and 100 GHz (from top 
to bottom). At each frequency, the solid line refers to the 
case of no-evolution, the dotted line to the pure luminos- 
ity evolution models for steep- and flat-spectrum sources 
described in Sect. 3.4.1 of Dunlop & Peacock (1990). 
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Fig. 7. Frequency dependence of the total emission (in 
terms of antenna temperature) from cluster sources for 
z = 0.15, the median redshift of our cluster sample. The 
dashed and dotted lines correspond to steep- and flat- 
spectrum sources, respectively, while the solid line shows 
the total. 



steep-spectrum (a > 0.5), so that the ratio of flat to steep 
spectrum sources is 0.35(+0.065, —0.063), consistent with 
the ratio (0.41) we found for a random sample of 1548 
field sources. The mean spectral indices are aflat = 0.13 
and a s teop = 0.87. 

The ratio of 4.85 GHz luminosity densities of flat- 
to steep-spectrum cluster sources is 0.38. Adopting the 
above mean spectral indices, such ratio becomes 0.15 at 
1.4 GHz. Using this ratio, we partition the total clus- 
ter radio luminosity as log L c i uste r,flat,i.4GHz — 22.48 and 
logL c i ust er,stocp,i.4GHz — 23.29 and extrapolate these lu- 
minosities to higher frequencies using the corresponding 
mean spectral indices. 

The big extrapolation in frequency required to reach 
the spectral range where most SZ surveys are carried out 
(y > 20 GHz) clearly makes our estimates, based on lim- 
ited information on the 1.4 to 4.85 GHz spectral indices, 
liable to large uncertainties. In fact, the low-frequency 
spectral indices may not apply up to much higher frequen- 
cies. Also, sources with strongly inverted spectra, such 
as extreme GHz Peaked Spectrum sources, undetected by 
low frequency surveys, may become important at high fre- 
quencies. An important test of our estimates is provided 
by the 30 GHz cluster survey by Cooray et al. (1998). 
These Authors have detected at > 4<r a total of 35 ra- 
dio sources, towards 56 clusters. If indeed, as argued by 
Cooray et al. (1998), most detected sources are physically 
associated to the clusters, we may derive their luminosities 
attributing to them the redshift of the associated cluster 



and estimate a mean 30 GHz luminosity function of cluster 
sources. Extrapolating to 30 GHz the mean surface den- 
sities of flat- and steep-spectrum sources in our control 
fields, we find that the fraction of background/foreground 
sources in their cluster fields is ~ 25%, consistent with the 
estimates by Cooray et al. (1998). The derived 30 GHz lu- 
minosity function was correspondingly scaled down by a 
factor of 1.25. In Fig. [5] this estimate is compared with 
the extrapolation to 30 GHz of our 1.4 GHz luminosity 
function, using the previously determined mean spectral 
indices for flat- and steep-spectrum sources. Although, due 
to the smallness of the 30 GHz sample, the uncertainties 
are large and the region of overlap is limited, the agree- 
ment is reassuringly good, supporting our choices for the 
mean spectral indices and indicating that we are not miss- 
ing important source populations. 

The contamination of the SZ effect by radio sources is 
more usefully expressed in terms of antenna temperature 



AT A (v, z) 



L 



flat 



L 



steep 



2^A 2 r c 2 lustor (l + z)3 



(3) 



where £fl a t(steep) — -^cluster, flat(steep), 1.4GHz (^(1 + z ) 

/1.4GHz)~ Qflat ( stoc p) , kb is the Boltzmann constant and 
f cluster is the physical radius of clusters, assumed here to 
be 1.73 Mpc. 

The results are shown as a function of the cluster red- 
shift, for 4 frequencies, in Fig. [B] where we have considered 
two cases: no-evolution, as suggested by the analysis in 
Fig. (solid lines), and pure luminosity evolution (Dunlop 
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& Peacock 1990; dotted lines). The latter are probably to 
be taken as upper limits, at least for the redshift range 
covered by our cluster sample (z < 0.4). The frequency 
dependence of the antenna temperature for z = 0.15, the 
median redshift of our cluster sample, is shown in Fig. [7| 

5. Conclusions 

A cross-correlation of FIRST radio sources with Abell 
clusters has demonstrated an excess surface density by a 
factor ~ 5 of radio sources within a projected distance, r p , 
of 0.1 Mpc (for h = 0.65) from the cluster center. At larger 
radii, the profile of the excess density can be described by 
a /3-model (Cavaliere & Fusco-Femiano 1976) with a core 
radius of ~ 0.70 Mpc and = 1.65 [see Eq. ©]. 

Within r p = 0.25 Mpc, chosen as a compromise be- 
tween the conflicting requirements of maximizing the 
sample size on one side and of minimizing the fore- 
ground/background contamination on the other side, we 
have a total of 617 sources within 951 Abell clusters, 
at redshifts $ 0.4, with a median surface density a fac- 
tor of 2.7 higher than in the field, implying that 63% 
of sources are cluster members. Their mean luminosity 
function does not show any evidence for cosmological evo- 
lution. Lumping all sources together and taking into ac- 
count the surface density profile we have obtained an esti- 
mate of the mean luminosity function at 1.4 GHz in excel- 
lent agreement with that derived by Reddy & Yun (2004) 
studying 7 nearby, X-ray luminous clusters, in the over- 
lapping luminosity range, and extending to luminosities 
two orders of magnitude higher. 

The lack of evidence for cosmological evolution up to 
z ~ 0.4 is at odds with the predictions of the pure lumi- 
nosity evolution models by Dunlop & Peacock (1990), in- 
dicating that either the evolutionary behaviour of sources 
in clusters and in the field differ or that such models over- 
estimate the radio source evolution at low redshifts. 

No significant differences between the properties of 
sources in clusters and in the field have emerged. The 
shape of our estimated luminosity function is very similar 
to the local luminosity function in the field at the same fre- 
quency (1.4 GHz), recently determined by Magliocchetti 
et al. (2002). The mean density is however about 3000 
times higher, a factor that corresponds to the matter over- 
density in clusters formed at z ~ 1.5. 

The distribution of spectral indices, derived exploiting 
the 4.85 GHz GB6 survey, is also consistent with being 
the same in clusters and in the field. Having estimated 
the fractions of steep- and flat-spectrum sources and the 
mean spectral indices of the two populations we have de- 
rived the mean emission of cluster sources, as a function of 
redshift and of frequency. The 30 GHz luminosity function 
obtained extrapolating our estimate at 1.4 GHz is found 
to be in good agreement with that derived directly from 
the 30 GHz observations by Cooray et al. (1998). 

If sources do not evolve, their mean contamination of 
cluster SZ effects, in terms of antenna temperature, de- 
creases with increasing redshift. For example, for obser- 



vations at 30 GHz with angular resolution matching the 
assumed cluster radius of 1.7 Mpc (h — 0.65) we esti- 
mate a mean contamination of ~ 13.5 /iK at z ~ and of 
~ 3.4 /iK at z ~ 0.5. The contamination increases by a fac- 
tor of 1.5 within 0.25 Mpc from the cluster center. If the 
pure luminosity evolution models by Dunlop & Peacock 
(1990) are adopted, the contamination level turns out to 
be very weakly redshift dependent. Obviously, the con- 
tamination level drops rapidly with increasing frequency 
and becomes pretty large at wavelengths above a few cen- 
timeters. 
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